Electronic Correlations in Manganites 
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The influence of local electronic correlations on the properties of colossal magnetoresistance mangan- 
ites is investigated. To this end, a ferromagnetic two-band Kondo lattice model is supplemented with 
the local Coulomb repulsion missing in this model, and is analyzed within dynamical mean-field the- 
ory. Results for the spectral function, optical conductivity, and the paramagnetic-to-ferromagnetic 
phase transition show that electronic correlations have drastic effects and may explain some exper- 
imental observations. 

PACS numbers; 71.10.Fd, 71.27.-fa 



The double exchange mechanism is generally con- 
sidered to be the origin of ferromagnetism in mangan- 
ites with perovskite structure, like Lai_a;Caa;Mn03 or 
Lai_2;Sra;Mn03. Due to the renewed interest in these 
compounds with respect to their colossal magnetoresis- 
tance (CMR) [||, the ferromagnetic Kondo lattice model 
(KLM, also known as s-d model) has recently been in- 
vestigated intensively |^|| as the microscopic basis of 
double exchange. Due to the crystal-field splitting of the 
five Mn d-orbitals into three degenerate t2g-orbitals and 
two degenerate Cg-orbitals at a higher energy level, three 
t2g-electrons can be approximately modeled as a spin- 
3/2, with the remaining n = 1 — x electrons occupying 
the eg-orbitals. The latter are delocalized via an effective 
double exchange transfer [Q t which is mediated by the 
overlap between Mn eg- and O p-orbitals. This trans- 
fer, together with the local Hund's rule exchange J > 
between t2g- and Cg-spin, constitute the ferromagnetic 
KLM 



KLM 



Here cj^^ and c^^^ are creation and annihilation operators 
for electrons on site i within Cg-orbital v with spin u; 

= kJ^aa' cta^aa'^zi^a' deuotes the Cg-spin (r: PauH 
matrices), the t2g-spin, and (ij) is the sum over all 
sites and their neighbors. At doping x, Yl,va'\^ivG^iva) = 
\ — X. To keep the model approach simple, the transfer t 
between electrons on neighboring Mn sites is assumed to 
be diagonal in the orbital index and independent of the 
lattice direction. 

For J ^ the eg-electrons are oriented parallel to the 
t2g-spins. With this constraint, a ferromagnetic align- 
ment of the t2g-spins allows for a maximal kinetic energy 
gain of the Cg-electrons. This double exchange mecha- 
nism is generally considered to be responsible for ferro- 
magnetism in manganites. However, it was pointed out 
by Millis et al. ^ that double exchange alone cannot de- 
scribe the resistivity of manganites. In particular, the 
importance of electron-lattice coupling was stressed [^[^] 
and the effect of lattice polarons and bipolarons was stud- 



C ■ C ■ 



2"^EE^-^- (1) 



=1 (y> 



ied |^,0j9| . The effect of Anderson localization , tun- 
neling magnetoresistance between phase separated do- 
mains and orbital polarons was also investigated. 
In spite of this, the unusual properties of the paramag- 
netic "insulating" phase of manganites which shows a 
(quasi-) gap in the photoemission spectrum and in 
the optical conductivity p^ ] for arbitrary doping x are 
not yet sufficiently understood. 

In this Letter we wish to emphasize the importance 
of electronic correlations due to the local Coulomb re- 
pulsion for understanding the properties of manganites. 
These correlations are neg lected in the KLM (|). This 
neglect is a priori unjustified in regard to the fact that 
the electronic repulsion is known to be the largest energy 
in the problem, and it originated only in the technical dif- 
ficulties of the ensueing many-body problem. Presently, 
the best manageable approach for this type of problem 
is the dynamical mean-field theory (DMFT) p4 15|. We 



use this approach here and show that (i) electronic cor- 
relations have a pronounced effect on the paramagnetic 
phase, yielding an upper Hubbard band in the spectral 
function and a shift of spectral weight as is observed ex- 
perimentally, and that (ii) the microscopic origin of fer- 
romagnetism changes from double exchange at doping 
0-5 < X < 1 to superexchange at a: — > 0. 

A correlated electron model taking into account the 
Coulomb repulsion between Cg-electrons has the form 
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Here, i^KLM is the KLM (|l|) for an Ising type t2g-spin 
of size \S^\ = 1 [|l6|, ni^cr = U and Vb are the 

on-site Coulomb repulsions within the same and different 
Bg-orbitals, respectively, and Fq is the Ising-component 
of the Hund's rule exchange coupling between Cg-spins. 
From band structure calculations JTzt a nd photoemission 
plus x-ray absorption experiments ]18| , P^ these parame- 
ters are estimated as C/, Fp ~ 8eV and Fq, J ^ leV. The 
bandstructure results for the e^-bandwidth W = 1 — 2eV 



1 



eV 112 



is somewhat smaller than that measured {W = 3 — 4 
Up to now, the Coulomb repulsion was taken 
into account only in two limits. Firstly, at a; = the 
correlated electron model maps onto an effective Kugel'- 
Khomskii-type model at strong Coulomb repulsion [ pO[ . 
Here, magnetism is obtained in second-order perturba- 
tion theory in t from superexchange. Secondly, the fully 
spin-polarized phase, where the correlated electron model 
reduces to a one-band Hubbard or t-J model in the or- 
bital degrees of freedom, was studied |^l| . Even the anal- 
ysis of the one-dimensional problem is very demanding. 

To solve Eq. (j^), we employ DMFT where the lattice 
model @j maps onto a self-consistent single-site problem 
for the Green function 
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with Q ^ = G (S: self energy), partition func- 

tion Z = ^5^^^-^ expy4[^'*,^' ,S'^g-l], 
Grassmann variables ^'j^^n ^'i^d ^'^^n Matsubara fre- 
quencies ujn — (2ri -I- l)7r//3, and 

Ar,^ ,s^g-'] v:.„ KL + '^JS'] ^..n 



shown in Fig. at n = 0.7, W = 2, U = = 6, and 
— 1. In units of eV these are realistic parameters 
for manganites if W is identified with its bandstructure 
value. For comparison the dotted line shows the results 
for the ferromagnetic KLM [U = Vq = Fq = {)) which was 
solved following Furukawa Q. Without local Coulomb 
repulsion, the spectral function consists of two bands 
corresponding to configurations with Cg-spin parallel and 
antiparallel to the t2g-spin. Some of these configurations 
contain doubly occupied sites. Thus a Hubbard band at 
the energy J7,Vb is expected to develop once this repulsion 
is taken into account. Indeed, the drastic changes shown 
in Fig.|l| illustrate this effect, which goes along with a 
shift of spectral weight from the lower band at the Fermi 
energy {uj — 0) to the upper band, for a; > (24| . Such a 
transfer of spectral weight is observed experimentally in 
the doped system |l^] . In view of the present results this 
effect may be attributed to electronic correlations. An- 
other effect of electronic correlations is the broadening 
of the bandwidth due to the imaginary part of the self- 
energy. This may account for the discrepancy between 
the bandwidth obtained by bandstructure calculations 
ipTf and the photoemission experiment ||l^. While elec- 
tronic correlations can explain some aspects of the man- 
ganite spectrum, the spectral weight at the Fermi level 
itself none the less remains large. We cannot exclude that 
the exact treatment of the three dimensional model may 
explain the unusual nature of the paramagnetic phase. 
But, at least within DMFT, electronic correlations alone 
can also not fully elucidate the unusual properties of the 
"insulating" paramagnetic phase of CMR manganites. 



The single-site problem is solved numerically by quan- 
tum Monte Carlo (QMC) simulations using a two band 
version of the Hirsch-Fye algorithm ||2^ . Employing the 
symmetry Sf — > —Sf,Ga G-a m the paramagnetic 
phase, only one of two local Ising configurations of (||) 
needs to be calculated. The mean value of spin-up and 
-down Green function for this Ising configuration yields 
the paramagnetic Green function. Self-consistency is ob- 
tained by iterating (|^) with the k-integrated Dyson equa- 
tion 
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where N^{e) = ^p^rj^^ {W /2Y - e'^ is the non- 
interacting density of states (DOS) with bandwidth W . 

The QMC algorithm yields the Green function at 
imaginary Matsubara frequencies. To obtain the Green 
function G,^cr{^) as a function of real frequencies an 
analytic continuation employing the maximum entropy 
method p3[ | is performed. The resulting spectral function 
Ai,cr{i^) — — ImG'i/cr(i^)/7r for the paramagnetic phase is 
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FIG. 1. Spectral function A{uj) {uj — corresponds to the 
Fermi energy) for the manganite model at doping x = 0.3, 
W = 2, U = 8, Fo = 1, Vo = U - 2Fo, and J = 3/2Fo. 
Electronic correlations yield a satellite band at about Vo and 
a shift of spectral weight from the lower band at the Fermi 
energy to higher energies. 

Next we calculate the optical conductivity. Within 
DMFT this quantity is obtained from the particle-hole 
diagram only, since, in this limit, vertex corrections to 
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the optical conductivity vanish 
employ the formula 



Following [|l5| we 



Rea{uj) = -- JdeN°{e) J duj' lniG,{uj')linGe{Lo' +uj) 
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as derived for the hypercubic lattice by Pruschke et al. 
|p5| . Here ImG^ is the imaginary part of 
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The result is shown in Fig. ^j. Already for [/ — V^ — Fq = 
(ferromagnetic KLM) the optical conductivity deviates 
considerably from the results for the non-interacting sys- 
tem (J = 0). This is due to the scattering of the Cg- 
electrons at the disordered t2g-spins which leads to a non- 
Fermi liquid paramagnetic phase with a finite imaginary 
part of the self-energy at the Fermi energy. Electronic 
correlations induced by the local Coulomb repulsion fur- 
ther reduce the optical conductivity at small frequencies. 
A second peak (absent for J = 0), which results from 
contributions to the particle-hole diagram with particles 
in the upper and holes in the lower band of the spectrum, 
is smeared out by the electronic correlations. 




KLM, T=0.125 

U=8, 1=0.2 

U=8,T=0.125 

J=0,T=0.125 



CO 



FIG. 2. Optical conductivity cr(tj) a± x = 0.3, W = 2, 
f/ = 8, -Fo = l, Vo = 6, and J = 3/2. Dotted line: r = 0.2; short 
dashed line: T = 0.125; solid line: KLM {U = Fo = Vo = 0); 
long dashed line: non-interacting system {J = U = Fo — Vo~0; 
line broadening 5 — 0.01). The optical conductivity of the 
interacting system is seen to differ considerably from that of 
a Fermi liquid. 

Finally, we investigate the instability of the para- 
magnetic phase against long-range ferromagnetic order. 
To this end, we calculate the magnetic susceptibility 
x(T) which diverges at the Curie temperature Tc as 
X{T) = a{T - Tc)-^ within DMFT. Fig. | shows the 
results as a function of doping x together with two ana- 
lytically tractable limits of DMFT: (i) double exchange 



(dashed line) as described by the ferromagnetic KLM and 
(ii) superexchange (cross) at a: = and strong coupling. 
The KLM [case (i)] neglects electronic correlations and, 
at J ^ t, leads to double exchange with an energy gain 
proportional to the kinetic energy gain of the Cg-electrons 
in a ferromagnetic environment, i.e., Tc (x t. At a; = 1 
the Cg-bands are empty and no kinetic energy is gained, 
i.e., Tc — 0. On the other hand, at a; = 0, the two 
spin-polarized Cg-bands are half filled, such that the ki- 
netic energy gain, and hence Tc, is maximal. In case 
(ii), DMFT corresponds to Weiss mean-field theory for 
the effective Kugel'-Khomskii-type Hamiltonian which, 
at a; = 0, predicts an instability against orbital ordering 
and an additional instability against ferromagnetic order 
at Tc = Zt^/{Vo - Fo) - Zt^/{Vo + 2J) (Z: number of 
nearest neighbors). Since t <^U,Vo the critical tempera- 
ture of the superexchange mechanism {Tc oc i^/Vb) is an 
order of magnitude smaller than that of double exchange 
(Tc cx t). 
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FIG. 3. Curie temperature Tc for the phase transi- 
tion from the paramagnetic (PM) to the ferromagnetic (FM) 
phase as a function of doping x. Dashed line: KLM with 
W = 2 and J — 3/2; squares: correlated electron model 
(^) which also takes into account the Coulomb interaction 
between eg-electrons {U = 8, Vb = 6, and Fo = 1); cross: 
Weiss mean-field theory for (^. The correlated electron 
model is seen to describe a crossover from double exchange at 
0-5 < a; < 1 to superexchange at a; — > 0. 

Solving the correlated electron model (^ for arbitrary 
X numerically within DMFT, a crossover from double ex- 
change to superexchange is clearly observed. At x > 0.5, 
the critical temperature is relatively well described by 
the ferromagnetic KLM, i.e., by double exchange. Here, 
double occupations are rare since there are only few Cg- 
electrons. With decreasing x (increasing number of Cg- 
electrons) the local Coulomb repulsion strongly reduces 
double occupancies and thus becomes more and more im- 
portant. The kinetic energy gain, and thereby double ex- 
change, is reduced and superexchange becomes effective 
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instead. This crossover from double exchange to superex- 
change yields a maximum in Tc in qualitative agreement 
with experiment. Note that without the coupling to the 
t2g-spin, i.e., in a two-band Hubbard model, no ferro- 
magnetism was observed for < n < 1 at values of Fo 
typical for manganites . 

In conclusion, our results show that electronic corre- 
lations are certainly important for understanding CMR 
manganites and cannot be neglected. In the paramag- 
netic phase, electronic correlations lead to the formation 
of an upper Hubbard band and to a shift of spectral 
weight into this band if the system is doped. This shift 
of spectral weight from a lower to an upper band is a 
genuine correlation effect and may explain a similar ex- 
perimental observation jl^] . Another aspect of electronic 
correlations is the broadening of the spectrum due to 
the imaginary part of the self-energy. This might be the 
reason why bandstructure calculations yield only about 
half of the spectral width compared to the photoemission 
experiment. The main result for the paramagnetic-to- 
ferromagnetic phase transition is that double exchange, 
as described by the ferromagnetic KLM, can explain fer- 
romagnetism in CMR manganites only for doping x > 
0.5. At lower values of x, the suppression of double oc- 
cupations by the local Coulomb repulsion becomes more 
and more important and leads to a crossover from double 
exchange to superexchange. This results in a maximum 
of the Curie temperature in qualitative agreement with 
experiment. 
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